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Abstract: We study some symmetry and integrability properties of four-dimensional Einstein- 
Maxwell gravity with nonvanishing cosmological constant in the presence of Killing vectors. 
First of all, we consider stationary spacetimes, which lead, after a timelike Kaluza-Klein re¬ 
duction followed by a dualization of the two vector fields, to a three-dimensional nonlinear 
sigma model coupled to gravity, whose target space is a noncompact version of CP^ with 
SU(2,1) isometry group. It is shown that the potential for the scalars, that arises from 
the cosmological constant in four dimensions, breaks three of the eight SU(2,1) symmetries, 
corresponding to the generalized Ehlers and the two Harrison transformations. This leaves a 
semidirect product of a one-dimensional Heisenberg group and a translation group as resid¬ 
ual symmetry. We show that, under the additional assumptions that the three-dimensional 
manifold is conformal to a product space M x S, and all fields depend only on the coordinate 
along M, the equations of motion are integrable. This generalizes the results of Leigh et al. in 
arXiv:1403.6511 to the case where also electromagnetic fields are present. In the second part 
of the paper we consider the purely gravitational spacetime admitting a second Killing vector 
that commutes with the timelike one. We write down the resulting two-dimensional action 
and discuss its symmetries. If the fields depend only on one of the two coordinates, the 
equations of motion are again integrable, and the solution turns out to be one constructed by 
Krasihski many years ago. 
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1 Introduction 

Exact solutions to Einstein’s field equations and to their supergravity generalizations have 
been playing an important role in many developments of general relativity, string theory and 
high energy physics [1]. For instance they can teach us a lot of insight into the theoretical 
aspects of general relativity such as the vacuum structure, uniqueness theorems and so on. 
The importance of exact solutions, though, is not limited to the classical situation, but extends 
also to quantum gravity. Indeed, much of our current knowledge on quantum effects in strong 
gravitational fields comes from the study of classical black hole solutions. 

Generically, the construction of exact solutions to general relativity is a notoriously dif¬ 
ficult problem, since the underlying field equations are a system of coupled nonlinear partial 
differential equations of second order. Nevertheless, one may hope that this system becomes 
integrable when some sufficient amount of symmetry is present. During late 1970s, a variety 
of independent groups have established the integrability of Einstein’s vacuum equations in 
stationary and axisymmetric systems. This includes the discovery of Backlund transforma¬ 
tions by Harrison [2] and by Neugebauer [3], and a Lax-pair representation by Belinsky and 
Zakharov (BZ) [4]^. In particular, the results of [4] have been generalized in many direc¬ 
tions, e.g. to the Einstein-Maxwell system [6, 7], five-dimensional general relativity [8] and 

'^Shortly before [4] appeared, Maison [5] was able to rewrite the stationary axisymmetric vacuum Einstein 
equations as a ‘linear eigenvalue problem in the spirit of Lax’, and noticed that this could ‘nourish some hope 
that a method similar to the inverse scattering method may be developed’. 
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five-dimensional minimal ungauged supergravity [9] (cf. also [10-12]). These techniques of 
integrable systems allow us to construct nontrivial new solutions starting from a given seed 
by adding solitons in a simple algebraic manner. Perhaps one of the most exciting recent 
achievements for the applications of these techniques is the inverse scattering construction 
of black objects admitting multiple horizons with various topologies in five-dimensional pure 
gravity (see e.g, [13] for a comprehensive review). 

Apart from the construction of numerous exact solutions, the underlying mathematical 
structure behind these integrable systems has been worked out by many authors. Prior to 
the studies of integrability, Geroch made a pioneering analysis of solution-generating methods 
for Ricci flat spaces in presence of a single [14] and two mutually-commuting Killing vectors 
[15] (see [16-19] for an electromagnetic generalization). In the presence of two commuting 
Killing vectors, the target space isometry does not commute with the internal SL(2,]R) sym¬ 
metry, giving rise to an inhnite affine Lie group called the Geroch group. Hauser and Ernst 
were able to prove the conjecture that any stationary and axisymmetric solution can be de¬ 
rived in principle from a Minkowski seed by the Geroch group [20]. Gosgrove addressed the 
interrelationships between several solitonic systems and the Geroch group [21]. Later on, 
Breitenlohner and Maison (BM) unraveled the group-theoretical structure of solitonic meth¬ 
ods from the standpoint of a Riemann-Hilbert problem [22]. A very nice general analysis of 
the relation between the BM group structure and the inverse scattering method of the BZ 
approach was recently given in [23]. In addition, a close connection to nonlinear sigma models 
has also been widely discussed [24]. 

In view of possible AdS/CFT (and many other) applications, one may thus ask whether 
similar integrability properties still hold in the presence of a cosmological constant, for ex¬ 
ample for stationary axisymmetric Einstein spaces in four dimensions. The introduction of 
a negative cosmological constant has a strong impact upon the spectrum of black holes. Of 
most prominent is that the horizon of an asymptotically AdS black hole can be a compact 
Riemann surface of any genus [25-28]. This is in contrast to black holes in asymptotically flat 
spacetimes [29, 30]. This is not the end of the story, since even more exotic possibilities exist, 
e.g. noncompact horizons with finite area [31, 32] (for generalizations to higher dimensions 
and further discussions of the physics of these solutions cf. also [33, 34]). One may thus expect 
a rich spectrum of black objects in presence of a cosmological constant, with many of them 
perhaps still to be discovered. It is clear that the integrability of stationary axisymmetric 
Einstein spaces would simplify enormously the construction of such solutions. A main ob¬ 
struction for this program is that the metric cannot be cast into the Weyl-Papapetrou form in 
the presence of a cosmological constant. It is therefore obvious that the techniques available 
in the absence of A cannot be straightforwardly applied. 

First steps in the investigation of the integrability properties with nonvanishing A were 
undertaken in [35-38]. These papers developed solution-generating techniques for the sta¬ 
tionary vacuum [36, 38] and electrovac [37] Einstein equations with a cosmological constant, 
in four [37, 38] and higher [36] dimensions. The cosmological constant leads to a potential in 
the dimensionally reduced system, breaking the symmetries of the original sigma model, and 
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thus the usual solution-generating techniques can’t be applied anymore. Still, some restricted 
formalism of solution-generating method is still applicable. In spite of this limited utility, 
they turn out indeed fruitful to generate some new exact solutions [37]. 

Here we shall make a first step towards a systematic investigation of the integrability 
properties of Einstein-Maxwell gravity with nonvanishing cosmological constant in four di¬ 
mensions, by extending the work of [38]. In the first part of this paper we consider stationary 
spacetimes which are described, after a timelike Kaluza-Klein reduction followed by a du- 
alization of the two vector fields, by a three-dimensional nonlinear sigma model coupled to 
gravity, whose target space admits an SU(2,1) isometry group. It is shown that the potential 
for the scalars, that arises from the cosmological constant in four dimensions, breaks three 
of the eight SU(2,1) symmetries, namely the generalized Ehlers and the two Harrison trans¬ 
formations. This leaves a semidirect product of a one-dimensional Heisenberg group and a 
translation group as residual symmetry. We show that, under the additional assumptions 
that the three-dimensional manifold is conformal to a product space M x S, and all helds de¬ 
pend only on the coordinate along M, the equations of motion are integrable. Subsequently, 
we consider the purely gravitational case and assume the existence of a second Killing vector 
that commutes with the timelike one, i.e., we focus on stationary and axisymmetric Einstein 
spaces. We write down the resulting two-dimensional action and discuss its symmetries. If the 
fields depend only on one of the two coordinates, the equations of motion are again integrable, 
and the solution turns out to be one constructed by Krasihski many years ago [39-41]. 

The remainder of this paper is organized as follows: In the next section, we discuss 
the integrability of the stationary Einstein-Maxwell-A system by assuming that the base 
space takes a product structure M x S and the target space variables depend only on a 
single coordinate. We derive the Reissner-Nordstrom-Taub-NUT-(A)dS metric by exploiting 
the Hamilton-Jacobi method. This generalizes the results of [38] to the case where also 
electromagnetic fields are present. In section 3, we address the integrability of the Einstein-A 
system, by assuming a second independent Killing vector, and derive an interesting class of 
solutions. Finally, we conclude in section 4 with some final remarks. An appendix provides 
an attempt of a higher-dimensional generalization. 

2 Einstein-Maxwell-A system 

In this paper, we focus on 3 -|- 1-dimensional Einstein-Maxwell-(A)dS gravity, with action^ 

^ {R - - 2A) , (2.1) 

and equations of motion 

- ^Rgf^u + = 2 {F^,aFJ - , VpF^'^ = 0 . ( 2 . 2 ) 

^We use the signature (—,+,+,+). The Ricci tensor is defined as + 

pp per _ pp per 
piz-*- ap per-*- up' 
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The Faraday tensor can be locally expressed in terms of a gauge potential as F = dA. 

We shall investigate the integrability properties of the stationary Einstein-Maxwell-A 
system, which is an extension of the work in [38]. 


2.1 Dimensional reduction 

Let us consider stationary spacetimes admitting a Killing field which is timelike at infinity. 
Applying the algorithm of Kaluza-Klein reduction along the timelike direction, the metic and 
the gauge field can be decomposed as 

ds^ = -e-'^(df + A:„dx")2 + ho.^dx'^dx^^ , A = B(dt + K„dx“) + 5„dx" , (2.3) 


where early greek indices refer to three dimensions, and the fields h^/S, Ka, Ba, 4> and B are 
t-independent. Here and in what follows, the indices a, 13,.. are raised and lowered by 
and its inverse. Then the effective three-dimensional Lagrangian derived from (2.1) becomes 


= Vh 


+ 2e^daBd‘^B 


+ Ka,pB){G^^ + K^^B) - 2Ae‘^ 


(2.4) 


where Gap = daBp — dpBa and Kap = daKp — dpKa- It is convenient to dualize the two 
vector fields to scalars, which can be implemented by adding to (2.4) a piece containing two 
Lagrange multipliers G and V' that ensure the Bianchi identities, 


£( 3 ) ^ £( 3 ) ^ 2Ge^^^daGp^ + (^ + GB)e^l^^daKp^ . 
Variation of (2.5) w.r.t. Kap and Gap yields 

= d^'ip + Gd^B - Bd-yG , 

and 


^a/3 ^ ^- . 

^/h ^ 


(2.5) 


( 2 . 6 ) 


(2.7) 


These equations express the field strengths in terms of the twist potential V' and the magnetic 
potential C. Plugging (2.6) and (2.7) back into (2.5) leads (after dropping a tilde on to 


£(3) = Vh - (Ja, J“) - 2Ae^ 


where we have introduced the notation 

1 


{Ja, jp) = - da^'dpcf) + de^’^UaU^p “ {daBdpB + daGdpG) 


( 2 . 8 ) 


(2.9) 


The equations of motion following from the Lagrangian (2.8) are the three-dimensional Ein¬ 
stein equations 

G^ll + ke^hap = {Ja, Jfi) - \hap (J^, J^) , (2.10) 
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supplemented by the divergence-type equations of motion 

Va[d^(p + 2e*{Bd^B + Cd^C) - = 2Ae‘^ , V«(e2'^a;“) = 0 , (2.11) 

Va(e'^5“S-2e^'^Ca;“) = 0, V„(e'^5“C + 2e^'^Sa;“) = 0. (2.12) 

(2.8) describes a nonlinear cr-model with pseudo-Riemannian target space coupled to Eu¬ 
clidean gravity in d = 3, with a potential. The latter breaks part of the target space isome¬ 
tries. 


2.2 Nonlinear cr-model and broken symmetries 

The target space $ of the scalars in (2.8) is a Bergmann space corresponding to a noncompact 
version of CP^ [24, 37, 42], namely it describes a coset space SU(2,1)/S(U(1,1) x U(l)), 
endowed with the metric 

dsl = gij{ip)d^Uif^ = d(j)‘^ + 4e2'^(d^ + CdB - BdCf - 4e^{dB^ + dC^) , (2.13) 

where (p^ = {(j), ip, B, C). One can easily verify that 

Rij = ’ Cjjkl = kl i D[Rjklm = 0 ■ (2-14) 

Here Rijkl and Cjjkl are the Riemann and Weyl tensors constructed from the target space 
metric Qij and the covariant derivative Dj. The Bergmann space is a special Kahler-Einstein 
manifold with negative curvature. The last equation of (2.14) is a differential characterization 
of a symmetric space, while the second equation implies a quaternionic structure [43]. 

The eight Killing vectors of <1> generating the isometry algebra su(2,1) are given by 


= i2 = Cd^ + dBj ^2, = -Bd^ + dc , 

^4 = —Cds + Bdc , Cb = —2d^ + ^'ipd^ + Bds + Cdc ■ 


^6 = + 
+ 


^{e-^ - {B^ + C^) f - 2^^ 


d. 




Ob - 

y2' 


C{e-^ - (R2 + c^)) - 2i>B 
^7 = -dBd^ + \2^B - C{e-^ - (R" + C")) 
+(e-<A + 3C^)dB + (4RC - 2^)dc , 
^8 = -ACd^ + [2V;C + B{e-^ - (r2 + c^))^ 
+ (4RC + 2^)dB + (e--^ + - 3B^)dc • 


R(e-'^ - (R" + C")) + 2V'C dc, 


d.i 




(2.15) 


The first hve Killing vectors represent infinitesimal transformations that are linear in the 
scalars and comprehend a twist transformation, two electromagnetic gauge transformations, 
an internal U(l) transformation and a scaling one. The remaining three are the most inter¬ 
esting, due to the nonlinearity in the fields, and they are usually called generalized Ehlers 
transformation (^g) [44] and two Harrison transformations (^ 7 , >^8) [45]. 
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In order to see that these Killing vectors indeed generate the SU(2,1) symmetry, let us 
define 

^2^ = + iis + *(6 - * 6 )] ! = -^[-(^7 + iis) + i{^3 - * 6 )], 

= ^{2C5 + i^i + 2iCe), Kll = i^l + K3^ E2^ = (2.16) 

771 3 _ ^ { TJ \ ZJ \ Z 7 2 _ / 771 1\* Z? 1 _ Z Z7 3\+ Z 7 2 _ / 771 3\* 

^3 — + '”21 , .C'l — —[^2 ) , ^3 — (,-C'l 1 , i^3 — \^2 ) , 

where Hi, H 2 are Cartan generators defined by 

Hi = "-ii - iCe , ^2 = ^(6 - 6^4 - 2 ^ 6 ), [i^i, ^2] = 0 . 

One can easily verify that these vectors Ei^ {i,j = 1 , 2 ,3) satisfy the su( 2 , 1 ) algebra 

W,Ek^] = 5k^Ei - 6/Ek^ , K/ = 0. (2.17) 

Note that the dependence of the scalar potential 

V{(j)) = -2Ae^ (2.18) 

on the dilaton (j) breaks the invariance under nonlinear isometries and scalings. It is easy to 
see that the latter is recovered if we admit a rescaling of A. The theory described by (2.8) is 
thus invariant only under SU(2,1)/Hi, where Hi C SU(2,1) is a subgroup generated by ^ 6 ) 
^ 7 , 1^8 corresponding to the Heisenberg algebra 

[?7, ^s] = 4^6 , [^7, Ce] = Ks) '^e] = 0. (2-19) 

The five unbroken generators close themselves to form another one-dimensional Heisenberg 
subalgebra in semidirect sum with 

[ 6 )^ 3 ] = - 2 ^ 1 , [ 6 >?i] = [?3,?i] = 0 , 

[^ 7 , ^4] = (a4)/6', [ 6 , es] = (ct5)/0 , ( 2 - 20 ) 

where i,j = 1,2,3 and 



0 

0 

0 

/2 0 0\ 

0-4 = 

0 0 1, <75 = 

0 10 



loo V 


The Heisenberg algebra (2.20) realizes the fact that the constant (p space constitutes a Nil 
manifold, viz, one can view the four-dimensional metric (2.13) as a Wick-rotated Bianchi-H 
universe. 

The well-known solution-generating techniques [23, 46, 47] based on group theory can thus 
not be applied in presence of a cosmological constant. Moreover, the broken symmetries are 
also a first sign of the loss of complete integrability, valid for A = 0 after another dimensional 
reduction [22, 48, 49]. This implies also the inapplicability of the inverse scattering method [4, 
23]. In what follows, we shall perform an analysis of some remaining integrability properties, 
extending the results of [38]. 
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2.3 Hamiltonian formalism and first integrals 

In the spacetime admitting a single Killing field, the sigma model still couples to the base 
space hap represented by three-dimensional Einstein gravity according to (2.10). Because of 
the intricacy of this system, we usually simplify the problem by assuming further symmetries. 
In the absence of A, the base space is decoupled from the sigma model by assuming an axial 
Killing field. More precisely, the metric without A can be cast into the Weyl-Papapetrou 
form, and the base space part can be obtained by quadrature once the sigma model on is 
solved. Unfortunately, this decoupling does not occur in the presence of A as we will see in 
section 3. 

In this section, we follow a different path to arrive at an integrable system. Along the 
lines of the argument in [38] , we consider the case in which the base space admits only a single 
degree of freedom. Now we suppose that h^p describes a warped product space M x S, with 
S a two-dimensional manifold. Moreover we assume that all the scalar fields depend only on 
the coordinate representing M. To capture this more conveniently, let us introduce another 
scalar field k that describes a rescaling of the three-dimensional metric h^p, 

hap = khap ■ ( 2 . 21 ) 

Absorbing the warp factor into k, hap can be taken to be an unwarped product, 

hapdx°‘dx^ = dcr^ -|- dH^ , (2.22) 


where dfl^ is the line element on S. Under these settings, every quantity depends only on 
a single valuable a. In this case the trace and the cnr-component of the Einstein equations 
(2.10) become respectively 


1 / d^k\ 

k Vdcr^y 




(2.23) 

(2.24) 


It is clear that the scalar curvature must be constant as a consequence of the fact that 
the r.h.s. of (2.23) depends only on a and the l.h.s. is independent of a. Without further 
resrictions we can thus take = 21 with I = 0,±1, so that S must be a maximally 
symmetric space, dfl^ = d9‘^ + f^{6)dip'^, where 


M^) = ^sin(\/I(9) 


/ 

sin 0 , ^ = 1, 

< 0, Z = 0, 

sinh 0 , I = —1 


(2.25) 


One obtains then a classical dynamical system with five degrees of freedom, with action 


S 


dak 2 


1 

2k‘^ Vdcry 


{Ja, Ja) +21- 2Ake‘^ 


(2.26) 
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For future convenience we introduce a new evolution parameter r defined by 

= dr. 


With the new potential V = 2K — ^ and oj = w,-, the action (2.26) can be expressed as 

S = f Ldr with a Lagrangian 


L = 


2 L 






(2.27) 


where a prime denotes a derivative w.r.t. r. It is easy to see that (2.23) is the constraint 
H = L + 2V = 0. It then turns out more convenient to pass to a Hamiltonian formulation 
rather than working in a Lagrangian description. After a Legendre transformation one gets 



e 





+ 


e / 2 2 

-^{pb+Pc 


‘^CpBP^ + 2BpcP^ + (H^ + 


+ F. 


(2.28) 


The solution of this dynamical system is highly linked to the existence of commuting constants 
of motion. The Killing vector fields of su(2,1) can be promoted to functions in phase space, 
realizing a Lie algebra isomorphism, by means of the substitutions^ 


^ Pipi: [t]'Ipb , (2.29) 

where {(^^} = {(p, "0, B, C} and i = 1,..., 8 . The minus sign in front of Ci reflects the fact 
that the infinitesimal generators and the corresponding charges obey the same algebra up to 
the sign of the structure constants^. The only nonvanishing Poisson brackets between the Ci 
and the Hamiltonian are given by 

{77,(75} =-277 +4A, {H,Ce} = dH^P - 8Aip , 

(77, Cr} = -dBH + SAB , {77, (Pg) = -477(7 + SAC . (2.30) 


Since the Ci do not depend explicitely on r, we find immediately that (7i, ( 72 , ( 73 , (74 are four 
constants of motion besides 77. Moreover if we define the modihed function C 5 = C^—dAr and 
use the constraint 77 = 0, we recover the constant of motion linked to a scale trasformation 

=-277 = 0. (2.31) 

dr 

^Our convention for the Poisson bracket is {A,B} = Q^^^dAiAdNB = ~ where 

P = i(T 2 is the symplectic form. 

'^This can be shown as follows. Let Qi = Qi{<l^,Pi) be first integrals obeying the Lie algebra {Qi,Qj} = 
f^ijQk and let us denote the corresponding Hamiltonian vector fields by OnQi- For any function 

F = F{q\pi) in phase space, we have a formula Vi^O mF = —{Qi,F}. It follows that for the vector field ViY = 

we obtain V^YOmF = -{{Q^,Qi},F} = -{Qi, {Q„ F}} + {Qj, {Qi, F}} = 
— [Vi, Vj]^dMF, where at the second equality we used the Jacobi identity. This establishes [Vi, Vj] = —f^ijVk, 
as desired. 










The modification of to is a consequence of the necessity to rescale also A in order to 
maintain invariance under scale transformations. 

The only nonvanishing Poisson brackets between the constants of motion read 

{C2,C^} = -2Ci, {C2,C'4} = C 3 , {C3,C4} = -C' 2 , (2.32) 

{C' 5 ,Ci} = -2Ci, {C5,C2} = -C' 2 , {C5,C3} = -C3. (2.33) 

Among Cl, ( 172 , 6 * 3 , 6 * 4 , C 5 and the operators composed of them, the maximal set of commuting 
first integrals is given by H, Ci, C 4 , Cf + C 3 , and we fix the values of these first integrals with 
four constants E, v, ATi, A' 2 , 

H = E, p^ = Av, Bpc - CpB = Ki , {pB + Cp^f + {pc - Bp^f = K 2 . (2.34) 

We want to use these equations to solve the system, so we shall set i? = 0 only at the end of 
the integration procedure. 


2.4 Integrability: RN-TN-(A)dS solution 

Using (2.34), the Hamiltonian can be rewritten as 

g ^ g ^ g 3^ g 20 ^ 

H = —pI - ^ (4p)" + W ^ 

and thus the electromagnetic and twist part has decoupled from the other fields. In order to 
solve the Hamilton-Jacobi equation 


H 



dk' d4>) 


dS 


we use the separation ansatz 


S = W{k,(t))-ET, 


(2.36) 

(2.37) 


which leads to 

e-<7 fdwV e-^ kdWy 
2 \ dfc y 2k‘^ \ d4> ) 8k'^ ^ 

(2.38) can be solved by defining the new variables x 
Charpit-Lagrange method. The result is 


{K 2 + 16vKi) + V = E. (2.38) 
= /se'^, y = e~^ and applying the 


W{x, y) = -^^/2ax — (^E{v^ + ax) — Gal — 12a‘^y^ 

+ ^(772 + 16uAri)arccot( ^ ^ ) , (2.39) 

ov \v2ax — v^J 

where E = 2A — E and a is an integration constant. Following the Hamilton-Jacobi technique 
we can introduce two other constants f3i, f32 according to 


/?! = 


55 

dE ’ 




55 

da 


(2.40) 
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Using the dynamical constraint = 0, they are given by 

,01 = 77 ^ \/‘lax - v^{v^ + ax) + r , 
oa^ 

^ A{2v^ — 2av^x — a?x^) ^ K 2 + Iv"^ — lax + 2a?‘xy 

'ia^y/2ax — 16aV2ax — a?y/2ax — v'^ 

To simplify the solution, it is convenient to define a new evolution parameter r by 

^ ( 


T = 


J.3 y2 
VT + 


(2.41) 

(2.42) 

(2.43) 


V 3 2 a^ 

To solve the two algebraic equations (2.41) and (2.42), we note that it is possible to set (3i = 0 
without loss of generality by shifting r. Then (2.41) gives 

2 

(2.44) 


2 

X = r + — . 

2a 


Plugging this into (2.42) yields 
1 


y = 


2 a(r^ + f^ 


Iv'^ A 


- v/2/32a3/2r + Ir^ -( r** + 

16 ^ 2a 3 


4 ^ 3r‘^v‘^ 3i;'^\ 


Using the original expression for H (2.28), the Hamilton equations for the electromagnetic 
part become 

dpc V 


dpB 

dr 

(W 

dr 




V 

+ g) 


{pc + AvB) 


1 


dr 


dC 


(r= + g) 


(-PB + 4oC), (2.45) 


4V2l (r= + g) 


{PB - dvC) , — = 


<!>■ (r= + g) 


{pc + dvB). (2.46) 


Using the gauge freedom generated by and ^ 2 ) we can implement a boundary condition in 
such a way that B and C vanish at infinity [42]. This eliminates two integration constants 
and the solutions are given by 


B = 


+ r/Sj 

2 , ’ 
^ + 2ii 


c = 


VTark-^^ 


_ 2q? 

r 2 + — 

' w 2a 


(2.47) 


Finally, the twist potential ^|J can be found by inverting the equation p^ = dv, which leads to 

dB „dC 


1/; = / dr — 




r^ + — 
' w 2a 




(2.48) 


The integration procedure is now complete. Since the constants defining the solution are not 
very illuminating, we define the new constants 


m = 


/32a3/2 

VT’ 


n = 




Q = 


2a = , (2.49) 


n 
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which give Ki = 0 and K 2 = 16(P^ + Q^)- It turns out that the four-dimensional metric and 
U(l) gauge field take the form of the RN-TN-(A)dS solution [50], 


ds^ = -e-^{dt + K^dipf + ke^ 


+ d9‘^ + {9)d(p'^'^ , 


(2.50) 


A^dx^ = Bdt + A^pdif , 


(2.51) 


where 


A = /(r^ — n^) — 2mr - {r^ + — 3n^) + , 

3 


k = 
B = 


k 


rrP + Pv? ’ 

Qr — nP 


e ^ = 


rz _|_ 


2 ’ 


Aip — 


= -4n^/m^ + Pn^ff {9/2 ), (2.52) 

2/2(0/2) (P(n2 - r2) - 2nQr) 


^/m? + PrP{r‘^ + rp) ’ rP + 

Note that the fields and are obtained from the dualization (2.6) and (2.7), that involves 

3lr — 3m — 4:ArPr\ ^ nQ + rP 


i) = 


n 


3(m2 -|- PrP) 


Ar -|- 


r2 -|- n2 


C = - 


y/rrpA/PprP{r‘^ -\- rP) 


For P = Q = 0 we recover the results of [38] , and thus the integrability properties described in 
[38] are still valid in the case of nonvanishing electromagnetic charges. We saw that, even if the 
cosmological constant reduces the internal symmetry group from SU(2,1) to SU(2,1)/Hi, it 
hasn’t spoiled integrability once we restrict to the subspace (2.22). This condition reduces the 
infinite number of degrees of freedom to effectively five. Only the three nonlinear generators 
of su(2,1) are broken and the remaining commuting first integrals are enough to decouple the 
electromagnetic and twist potentials and to integrate the system in three steps. The general 
case remains unsolved and is highly linked to the broken affine Kac-Moody algebra arising 
after another dimensional reduction [49, 51, 52]. The action of SU(2,l)/Hi on the fields 
generates a transformation on the parameter space, and in particular a scale transformation 
requires a rescaling also of A. Unfortunately these surviving symmetries alone are useless to 
produce new interesting solutions. 


3 Einstein-A system 

In this section, we shall consider the action (2.1) with vanishing electromagnetic field = 0, 
and assume the existence of an additional Killing vector that commutes with dt- For A = 0, 
this system is described by Weyl-Papapetrou formalism which allows us to utilize certain 
integrability techniques. We see that A term destroys the reduction to the Weyl-Papapetrou 
system. In spite of this, a further reduction to d = 1 with a suitable choice of variables enables 
us to solve the Einstein-A system in full generality. 
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3.1 Effective field theory in two dimensions 


In general, a solution with an M x SO(2) isometry group cannot be written in the Lewis form 
[53], but this becomes true if the line element admits a two-dimensional foliation orthogonal 
to the one in which the action of M x SO(2) is transitive. With this additional hypothesis we 
take 


ds^ = —e '^(dt -|- Kdip)'^ + e'^(e^'^Lmr).dx'"dx"' -|- , 


where t and (p are Killing coordinates, and the metric depends only on x'” {m 
Plugging this into (2.1) (with = 0) yields the two-dimensional Lagrangian 


£ 0 ) _ y^gX 




(3.1) 
1 , 2 ). 

(3.2) 


12 ( 2 ) 

is the two-dimensional Ricci scalar associated with hmn- In order to emphasize the 
nonlinear SL(2, M) symmetry present for A = 0, we define new fields il^, by 2^ = —cj)+ax+‘^ip 
(where a is an arbitrary constant to be specified later), (p = (j) — x, which leads to 


£( 2 ) = ^gX 


12 ( 2 ) + 2dkXd’"^ + 



dkXd^X 


(Jfc, J^) - 2 Ae"^+ 2 ^ 


(3.3) 


where the cr-model target space is SL(2,M)/SO(2), and thus 

{Jk, j^) = lidkPd'^P - e-^hiKd^K). (3.4) 

When A is turned off, we can recover another SL(2,M) via dualization in d = 3 [22, 24, 42]. 
Since these two nonlinear representations of SL(2,M)^ are linked by a nonlinear and nonlocal 
relation, they do not commute. This non-commutativity can be used to generate new solutions 
[46, 47] and to study the symmetries of the effective field theory (3.3) (with A = 0) [22, 48]. 

Varying (3.3) and using the fact that in two dimensions every metric is conformally flat®, 
one obtains the equations of motion 

Ae^ + 2Ae(i+")>^+2^ ^ g , (3.5) 

diie^d^p - e-^^+^Kd^K) = 0 , K) = 0, (3.6) 



Jz) - {Jq, Jq)) - d^xdz^ + dgxdg^ 
-^ (^ + a) {dzXdzX - dgxdgx) - Ae"^+2b 


(3.7) 


dzdgC^ = 


dzXdg'P + dgxdz'ip + 



dgXdzX 


{Jg, Jz) ) 


(3.8) 


®The first SL(2,R) is called the Matzner-Misner [54] transformation and the second one the Ehlers trans¬ 
formation [44]. 

®The conformal factor can be absorbed into p), cf. (3.1). 
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+ “) + 2 = 0 > (3-9) 

where A = + 3^. It is interesting to note that for the choice a = 0 (3.5) becomes the 
massive Klein-Gordon equation on a curved space with metric e^^(dp^ + dz^), and mass 
m? = —2A. In the case A = 0 it boils down to Ae^ = 0, which allows to take = p without 
loss of generality, but for A 7 ^ 0 y does not decouple anymore from the other fields. This 
explains the failure of the metric to fall into the Weyl-Papapetrou class. Notice also that 
(3.9) follows from the other equations. To see this, consider dz{3.7) - 9^(3. 8 ) and use (3.5) 
and (3.6). This leads to 

^d^xidlx - d^x) + - dzi’dkxd'^X = 

a dkxd^xdzX - {Jz, J^)dkX - dzXi^^ “ Aae“^+^’^). (3.10) 

Moreover one can rewrite (3.7) - {3.8)-dpX in the form 

^dzxidzX - dpX) + ^dzidgxf - dzi’dkxd'^x = 

(f ~ i) - {Jg, Jz)dgX + \dzX{{Jg, Jg) “ {Jz, Jz)) ■ (3.11) 

Finally, the difference of (3.10) and (3.11) implies (3.9). 

Note that the cr-model isometry group SL(2, M) acts only on cj) and K. Since the potential 
in (3.3) is independent of 4> and K, this SL(2,]R) is a symmetry of the complete Lagrangian 
(3.3). However, one easily checks that, from a four-dimensional point of view, a transformation 
with this SL(2,M) corresponds merely to a diffeomorphism, and can thus not be used to 
generate new solutions. 


3.2 Further reduction to d = 1 


If we make the additional assumption that all the fields depend only on g, the two-dimensional 
effective field theory (3.3) boils down to a dynamical system with four degrees of freedom 
described by 


S 


dge^ 


„dY d^ 




1 ( # Y /di^V 

2 Vd^iJ 


2Ae“x+2b 


(3.12) 


which turns out to be exactly solvable. Introducing the new coordinate r by e ^dg = dr and 
defining V = 2Ae^'^, the system (3.12) is described by the Lagrangian S = f Ldr as 

L = 2x'^' + (^ + «) , ( 3 . 13 ) 

where a prime denotes a derivative w.r.t. r. In what follows, we shall make the choice a = —2, 
for which x becomes cyclic. Then the equations of motion following from (3.13) are given by 

^// _ ^ Q ^ ^ g20^ ^ 
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= 0 , x" + 2Ae2^ = 0 , (3.15) 

where A is an integration constant, together with the constraint H = L + 2V = 0, that 
emerges from (3.7). Eqs. (3.14) and (3.15) are easily solved', and one finds that the metric 
in four dimensions is given by 

ds^ = + Kdcp)^ + e^^dr^ + + e'^~^^d(p^ , (3.16) 


where 


e '^ = -^= cosl^/^(r + C 2 )] 

grV4C3+3Ci/3 
^ [cosh(-v/C3r)]2/3 ’ 



C 3 

3Acosh^(\/C3r) ’ 
tan[v^(r- + ^ 2 )], 


(3.17) 


and Cl, C 2 , C 3 , C 4 , are integration constants. For generic values of these constants, dt, 
and dz constitute an exhaustive list of Killing vectors, as one can verify by checking the inte- 
grability of the Killing equation. This solution falls into Petrov type I and the Kretschmann 
invariant is 


= ^^/2 [3(^3cosh^( v^?-) + Cg^^ sinh'^(yC^r) 

3 C 3 

-2(C3 + 3Ci) v'4 C3 + 3Ci cosh3(v^r) sinh(V^r)] , 


while the Chern-Pontryagin one vanishes. Hence r —)> 00 corresponds to a curvature singu¬ 
larity. This class of solutions has already been found by Santos and MacCallum in [55, 56] 
and previously by Krasihski in [39-41], in both cases in a different coordinate system. The 
connection between the metric (3.16) and those of [55, 56] is easy to find after the iden¬ 
tification = 'k(r). Other explicit forms of the fields defining (3.16) can be recovered 
through the combination of analytic continuations, diffeomorphisms and nontrivial limits on 
the constants, requiring the final metric to be real. 

An example of a well-known static solution inside the class of (3.16) is found if we define 
A = \/Cia and make the choice 

Cl = 0 , 0 = 1 , C 3 = 9M ^, C 4 = , (3.18) 

where the parameter £ is related to the cosmological constant by A = —3£“^. Then the line 
element boils down to 

ds^ = -e^dt^ + e^d(^2 ^ ^ ^2{i^-x)^^2 ^ 


with 


e 2 b = 


M 2£2 


cosh2(3Mr) 


= (M^2)2/3. 


^2Mr 


cosh(3Mr)]2/3 


(3.20) 


^The case A — Q leads, after a change of coordinates, to (22.8) of [1], and it cannot be recovered smoothly 
after the integration. 
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Introducing the new radial coordinate R by 


R = (M£2)V3. 


^Mr 


cosh(3Mr)]i/3 ’ 


(3.21) 


the solution becomes 


ds^ = R^{-dt^ + ^ 1 dR^ + ( ^ 1 dz^ . 


2M . R^\ 


R £2 


2M 


R £2 


(3.22) 


This is the well-known planar AdS soliton which turns into a planar black hole after the 
analytic continuation t ^ iz, z ^ it. 

Another interesting limit of (3.16) is obtained if we choose C 3 = 3A, Ci = (3“^, and, after 
the limit A — )• 0 and the rescaling (r, z) 1—)> (-v/C^r, -^/^z), we fix the other parameters as 

A = /3 = yj, C2 = 0, C4 = ^. (3.23) 

Then the metric (3.16) becomes 


k^ds^ = —e'’[cos(\/3r)(dt^ — dp^) + 2 sin(\/3r)d(^dt] + e ^’"dz^ + dr^ . (3.24) 


This is the Petrov solution [53], eq. (12.14), also found as the metric induced on a constant 
r hypersurface of the d = 5 Einstein-Maxwell-A solution discussed in [36], eq. (5.37). It 
is interesting to note [1] that (3.24) is the only vacuum solution of Einstein’s equations 
admitting a simply-transitive four-dimensional maximal group of motions generated by the 
Killing vectors [53] 

T = dt, Z = dz, ^ = d^, R = dr + zdz + - (f) + V^p) dt , (3.25) 

satisfying 

[R,T] = ^T-^^, [R,<^] = h + ^T, [R,Z] = -Z. (3.26) 

Furthermore the determinant of (3.24) is always —1, so {t,4>,r,z) can be promoted to global 
coordinates. 

It is worthwhile to note that in our approach, which stresses the symmetries of the 
underlying theory, a simple system of integrable equations (3.14), (3.15) emerges in a natural 
way. However finding the most general solution to (3.5)-(3.8) remains a hard problem, due to 
the breaking of the SL(2,M) Ehlers symmetry, similar to the Einstein-Maxwell-A case. As we 
tried to argue in the introduction, this system of equations may nevertheless admit a Lax-pair 
representation, which we expect (if it exists at all) to be highly nontrivial to find. We hope 
to come back to this point in a future publication. 
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4 Final remarks 


The integrable nature of Einstein’s gravity in certain contexts is a useful ingredient for our 
understanding of the nonlinear nature of gravitational physics. In this paper, we developed 
some novel techniques of integrability that allowed us to obtain solutions of the Einstein- 
Maxwell-A system. In the first part we derived the complete solution by assuming that every 
quantity depends only on a single variable and that the specific form of the base space is M x , 
by extending the work of [38]. This restriction is strong enough to give a tractable system, 
yet turns out rich enough to incorporate a class of gravitational solutions of physical interest. 
In the second part, we found that under the codimension one assumption, an appropriate 
choice of variables allows us to obtain the decoupled equations (3.14), (3.15). We expect that 
these convenient variables would continue to be powerful for further investigations of the case 
with higher codimension. 

The present work can be generalized into various directions. An interesting plausible 
route is to see if other gravitational theories display integrability properties similar to the 
Einstein-Maxwell-A system in four dimensions. As we shall discuss in appendix A, it seems 
that the higher-dimensional generalization is not straightforward even for the pure Einstein- 
A case. Nevertheless, this failed attempt also gives further insight into the integrability 
properties of Einstein’s equations. 

Another possible future work is to extend our formalism to the case where the base space 
admits the dependence on more than one variable. To this aim, hints that these systems could 
actually be integrable come from the study of the Plebahski-Demiahski solution [57], which 
is the most general known Petrov-type D solution of the four-dimensional Einstein-Maxwell 
equations with A. It possesses a very high degree of symmetry and contains a lot of subcases 
and parameters of physical interest, so one may think of being able to generate it from a given 
seed solution. Also its relation to supersymmetry [58] may be a tractable way to understand 
the integrable nature of Einstein’s gravity in the presence of a cosmological constant. We 
hope to come back to these points in future work. 
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A On the integrability of the higher-dimensional Einstein-A system 

In this appendix we discuss if the integrability properties in four dimensions can be extended 
to higher dimensions. 

Let us consider the D-dimensional spacetime {Mo^go) represented by Einstein’s gravity 
with a cosmological constant, 

Ad = ^ / - 2 A) . (A.l) 
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Suppose that {MD,gu) admits n Killing vectors which mutually commute. The metric can 
be put into the form 

dsl) = gij.u{x)dx>^dx'' + Mmn{x){dy^ + + 2K^”^(rE)da:''). (A. 2 ) 

Here 5 ^ 1 / is the metric of the external d = D — n dimensional space, M^n is the internal 
metric and are the Kaluza-Klein gauge fields. Due to the isometries, we can reduce 

the D-dimensional gravity system down to d dimensions. Denoting the Kaluza-Klein field 
strength as it is straightforward to show that the Ricci scalar is 

decomposed as 

RiD) =R _ 

- , (A.3) 

where is the linear connection compatible with g^y. Let us perform the conformal trans¬ 
formation 

9i,y = ke‘^‘^g^y , Mrnn = Mmn , det(M) = -s, (A.4) 


where s = ±1 and the additional scalar field k has been introduced as in the body of text. 
To achieve the Einstein frame for A: = 1, we choose 


2 n 


a = 


{d + n-2){d-2) ’ 


/? = - 


{d — 2 )q; 


n 


(A.5) 


One can then verify that the d-dimensional equations of motion can be derived from the 
action 






R - + iTr(VMVM-i) + \{d- l){d - 2 )^^ 




- 2Afce“‘^ 


(A. 6 ) 


Hence M^n is the matrix corresponding to an SL(n,M)/SO(n) nonlinear sigma model. 

In the following, we shall focus on the d = 3 case. This is a generalization of the vacuum 
[59] and electrovac [60] analyses. The Maxwell equations can be solved to give local twist 
potentials V’(m) such that 

, (A.7) 

where e^yp is the volume element compatible with the external metric g^y. It therefore turns 
out that the 3-dimensional Einstein- and matter field equations can be derived from the action 


‘S's = / d^x^/\g\k^/‘^ 


2^3 


R -h |Tr(VMVM“^) -F i^i'^kf - 2AA:e“'^ 


(A. 8 ) 
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where 


(detM) ^ — (detM) 

-(detM)"V(n) Mmn + (detM)“ 


, detM = 1, M = , (A.9) 


with inverse 


.1 ^ / detM + MP'?V’(p)V’(p) M-PV'(p) \ 
M”^V'(p) M”"” y 


(A.IO) 


It turns out that the coset symmetry is now enhanced to SL(Z) — 2)/SO(L) — 2). 

Let us focus on the D = 5 case in what follows. Assuming two commuting Killing vectors 
^( 1 ) = d/dt and ^( 2 ) = djdip, we parametrize 


Mmudy^dy"^ = -f{dt + u^dif^f + (fh) , 
It follows that the matrix M in (A.9) can be expressed by 

/iV’(i) 




( -h /iV’(i) /iV'(2) 

^V’(i) -/ - -/w - hV'{i)V’(2) 

\ ^V’(2) - H{i)'4^{2) -HI2) - 


(A.ll) 


(A.12) 


This coincides with the parametrization {qi,q 2 ,Q 3 ,Pi,P 2 ,P 3 ) of SL(3,M)/SO(3) given in [61] 
by the identification 


qi^-h, q2^-f, Pi -'0(1) , P2^u), P3'-^-0(2)- (A.13) 

The Killing vectors Xi {i = 1, ..., 8 ) of the target space ds| = —(l/4)Tr(dMdM“^) are given 
by (5.24) of [61]®. 

For k = 1, the potential depends only on 0. The Killing vectors which keeps the potential 
invariant are given by 

X_ = Xi- 2A2, A3, A4, A5, Ag. (A. 14 ) 

The remaining three generators (A_|_ = A 2 — Ai, A 7 , Ag) do not leave the potential invariant, 
and form a Sim(l) x M algebra corresponding to the Bianchi III universe, 

[A+, A7] = A7 , [A+, As] = [A7, As] = 0 . (A. 15 ) 

We make a codimension one ansatz and require that the base space metric has an SO(3) 
symmetry. Then, the base space reads 

ds| = dx^ + 0-4 + (T 2 , (A. 16) 

seems to have a typo in ref. [61] and it should be modihed to Xq = 2 q 2 P 2 dq 2 d-psdpi + —pV)dp 2 ■ 
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where ai are SU(2) invariant forms, 


(Ti = — sin ^jJdO + cos V’ sin Odcp , a 2 = cos 'ipdO + sin V’ sin Od^p , us = d'll) + cos Odp . 


Every quantity is dependent only on a single variable x. This class of metrics includes the 
five-dimensional Myers-Perry black hole [62] with equal angular momenta. Defining 

k^/^hdx = dr, V = 2K- 2k-^h -^, (A.17) 


the one-dimensional Lagrangian S = f Ldr reads 

L = ^hk'^ + jA:2/iTr[(M)'(M-^)'] - V, 

where the prime denotes differentiation w.r.t. r. The Hamiltonian boils down to 
1 


(A.18) 


H =- 
2 L 


h Vi + / '^k '^p'i + f 'h '■^k '^pip^-fh '^k {Pi,, + uip^^) 


0 - 2 , - 2,-2 2 


-U-3,-2 2 


,- 2 ;- 2 / 


^p)- fpfPh + hpl) 


+ E. 


(A.19) 


The trace of Einstein’s equations gives the constraint H = d. For = {/, w, h,' 0 ( 1 ), V’( 2 )}) 
one can verify that are obvious first integrals and other nonvanishing Poisson 

brackets with the Hamiltonian are 




{C2,H} = -‘^H+^K, 


{Cr, H} = -20(1)77 + 4A0(i), {Cg, H} = -20(2)77 + 4A0(2). 


(A.20) 

(A.21) 


Hence C_ = Ci — 2 C 2 is another first integral. The nonvanishing Poisson bracket among 
these first integrals read 

{C'3,C4} = -C5, {774,C6}=C'_, {C5,C<i} = C3, {C_,C3} = C3, 

{C'_,C4} = -2C4, {C_,C5} = -C5, {C.,Cq} = 2C6. (A.22) 


Thus C 4 ^,Cq,C- generate an SL(2,M) subalgebra and its quadratic Casimir is 

another obvious first integral. Clearly, this is not true for the SL(3,M) quadratic Casimir 
\{Cl — C 1 C 2 + C 2 — C 4 CQ — — C^Cg) which involves also generators of the broken 

symmetries. 

We found that three is the maximal number of commuting first integrals in involution 
with 77 and there are many sets of this type that can be built from C-, C 3 , ( 74 , 6 * 5 , Cg and 
their compositions. In order to decouple some of the fields in the Hamilton-Jacobi equation 
associated to (A.19), the most intriguing set seems to be composed of 77 ,( 73,(75 and the 
cubic invariant = ( 7|(74 -I- C^Cq — CsC^C-. Unfortunately, we were unable to find a 
sufficient number of commuting first integrals to carry out the integration procedure as in the 
four-dimensional Einstein-Maxwell system. One plausible reason for this is that the target 
space SL(3,M)/SO(3) is five-dimensional admitting eight Killing vectors, while in the 77 = 4 
Einstein-Maxwell case the target space (2.13) is four-dimensional with eight Killing vectors. 
Obviously, the former is less symmetric. It would be interesting to see if this system displays 
chaotic behavior. We shall leave this point for future investigation. 
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